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ABSTRACT 

Background: Bovine theileriosis results from infection with obligate intracellular protozoa of the 
genus Thai ffia. The phylogenetic relationships between two isolates of Thdleria annulata, and 36 Thal& 
ria spp., as well as 6 outgroup including Babesia spp. and coccidian protozoa were analyzed using the 
18S rRNA gene sequence. 

Methods: The target DNA segment was amplified by PGR. The PGR product was used for direct 
sequencing. The length of the 18S rRNA gene of all Thdleria spp. involved in this study was around 
1,400 bp. 

Results: A phylogenetic tree was inferred based on the 18S rRNA gene sequence of the Iran and 
Iraq isolates, and other species of Thdleria available in G enBank. In the constructed tree, Thdleria an- 
nulata (Iran vaccine strain) was closely related to other T. annulata from Europe, Asia, as well as T. 
lestoquardi, T . parva and T . taurotragi aU in one clade. 

Conclusion: Phylogenetic analyses based on small subunit ribosomal RNA gene suggested that the 
percent identity of the sequence of Iran vaccine strain was completely the same as Iraq sequence 
(100% identical), but the similarity of Iran vaccine strain with other T . annulata reported from China, 
Spain and Italy determined the 97.9 to 99.9% identity. 
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Introduction 

Bovine theileriosis results from infection 
with obligate intracellular protozoa of 
the genus Theileria. The two most im- 
portant species in cattle and water buffalo are 
T. parva, as the agent for causes east coast fe- 
ver, and T. annulata, which causes tropical 
theileriosis. Theilffia parva occurs in 13 coun- 
tries in sub-Saharan Africa causing East Coast 
fever (ECF), Corridor disease, and January 
disease. Theilffia annulata, the cause of tropical 
theileriosis, occurs in large parts of the Medi- 
terranean coast of north Africa, extending to 
northern Sudan, and southern Europe. South- 
eastern Europe, the near and Middle East, In- 
dia, China and central Asia are also affected 
(1-3). 

Thalffia is classified based on microscopic ob- 
servations, ultrastructural features, lifecycle, 
geographic region and vertebrate and non- 
vertebrate host (4). Clinical classification cate- 
gorizes Theilffia to malignant, moderate and 
benign species. Nonetheless, the exact taxo- 
nomic Theilffia spp. have been difficult to es- 
tablish and the subject of considerable debate 
(5-7). There is a number of factors makes 
complexity of assigning taxonomic positions, 
including similar morphology among this 
group of parasites regardless of vertebrate 
host, incomplete lifecycle data, serologic tests 
that are not specific enough to discriminate 
individual species in the presence of mixed 
infections, and the difficulty of obtaining pure 
isolates for studies when the circulating para- 
sitemia may be very low (8-10). 
The advances in molecular biology and spe- 
cially sequence data analysis allowed the re- 
searchers to identify and characterize the he- 
moparasites species in particular Thalffia 
group. Ribosomal RNA is the most abundant 
constituent of nucleic acids in any non-viral 
organism with the eukaryotic RNA transcrip- 
tion unit consisting of the large and small sub- 
unit (18S rRNA) and the 5.8S rRNA gene (11). 
The 18S rRNA gene is increasingly accepted 
as a widely used marker for characterization. 



taxonomic classification, and phylogenetic 
analysis and this gene has been sequenced 
from a variety of different organisms, resulting 
in a large database for sequence comparisons 
(6, 10, 12-22). The conserved function and 
structure of the 18S rRNA molecule allow se- 
quences to be aligned, even among divergent 
species. However, the molecule also possesses 
phylogenetically informative variable regions 
that are useful for determining relationships 
among species (23). 

In this study, we decided to determine the 
phylogenetic position of T. annulata infected 
cell hne S15 Iran vaccine strain and T . annulata 
Iraq field isolate with 36 different Thaleria spp. 
and six out-group protozoan parasites using 
18S rRNA gene sequences comparison. 

Materials and Methods 

Paiasites 

Two T. annulata strain/ isolates were used in 
this study. Theilffia annulata infected cell Hne 
(S15 Iran vaccine strain, Tehran, Vasfenard) 
and Iraq T. annulata field isolate were provided 
from Protozoology and Vaccine Production 
Department of Razi Institute and Duhok 
Province, Kurdistan of Iraq respectively. 

DNA isolation 

Protdnase K and furtiier phenol chloroform 
purification were performed for DNA extrac- 
tion (24). Briefly, after treating the cells with ly- 
ses buffer, followed by centrifugation, pratei- 
nase K and SD S solution was added to the pellet, 
and then was incubated until most of the cellular 
protein was degraded. The digest was depro- 
teinized by phenol/ chloroform/ isoamyl alcohol 
extraction, recovered by etiianol precipitation, 
then was dried and resolved in TE buffer D NA 
concentration was determined either by agarose 
gel electrophoresis and spechDphotometry 
(A260) and measuring the ratio of A260/A280. 
Moreover, quality of the isolated DNA was eva- 
luated by agarase gel electrophoresis. 
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PC R Primer D esign 

The specific primers were designed based on 
T. annulata 18S rLbosomal RNA gene sequence 
(accession # EU083801) (by CinnaGen, Iran). 
Two primer pains were designed in order to 
span the major hyper variable regions along 
the IBS ribosomal RNA gene sequence (Gene 
Runner program, Version 3.05). The first two 
primers, Fl (5' GGC GGC GTT TAT TAG 
ACC 3') and Rl (5' TCA ATT GOT TTA 
AGT TTC AGC C 3') were used to amplify 
bases between 186-1093 and tiie second pri- 
mers, F2 (5'CAG ATA CCG TCG TAG TCC 
3') and R2 (5' CCT TGT TAG GAG TTC 
TGG 3') were applied to amplify bases be- 
tween 945-1714 of T. annulata IBS ribosomal 
RNA gene sequence (EU083801) and these 
two primers sets covered the majority length 
of IBS ribosomal RNA gene sequence with 
127 bp overlapping. 

Polymerase Chain Reaction 

PGR was performed in a final reaction volume 

of 20 pil containing 50 mM KGl, 10 mM Tris- 
HGl (pH B.3), 1.5 mM MgGl^, 0.1% Triton X- 

100, 200 [xM (each) deoxynucleoside triphos- 
phate, 0.5 U of Taq DNA polymerase (Ginna- 
G en, Iran), 10 pmol of each primers, and 2 pil 
of DNA template. The reactions were per- 
formed in an automatic DNA thermal cycler 
(Techne, G ermany) with the first incubation at 

94 G for 3 min and were pursued by 35 cy- 
cles. Each cycle consisted of a denaturing step 

of 20 seconds at 95 G, an annealing step of 45 

seconds at 56 G, and an extension step of 50 

seconds at 72 G, followed by final extension 

step of 10 min at 72 "c. 

PCR product detection and sequencing 

Amplified PGR products were separated by 
electrophoresis on a 2% agarose gel, stained 
with elliidium bromide, and visualized by UV 
transillumination. PGR products were cleaned 
and extracted from agarose gel and were sub- 
mitted for bidirectional DNA sequencing by 
using chain termination method (MWG , G er- 



many). The provided sequences from FlRl 
and F2R2, first were merged and offered to be 
aligned for multiple sequence alignment and 
phylogenetic study of Thdleria spp. and out- 
groups as well. 

Sequence alignment and phylogenetic 
analysis 

The DNA sequences of IBS rRNA gene ob- 
tained from two studied T. annulata samples 
and 42 sequences of IBS rRNA gene se- 
quences including 36 Thaleria spp., 4 Babesia 
spp. and two coccidian protozoa were ac- 
cessed from GenBank. The sequences were 
ahgned by Glustal W multiple alignments pro- 
gram (25). The alignment was manually edited 
in Bio Edit and truncated to the size of the 
smallest sequence (~l,400bp). Phylogenetic 
tree was constructed by using DNADist 
Neighbor-Joining method (version 3.6a2.1), 
sequence identity matrix of all sequences were 
computed as well (BioE dit phylogeny package. 
Version 7.0.1). 

Nucleotide sequence accession numbers 

The IBS rRNA gene sequences of the T. annu- 
lata S15 Iran vaccine strain and T. annulata 
field isolate from Iraq, Duhok have been sub- 
mitted to G enBank and can be retrieved under 
accession numbers of HM628581, and 
HM62B5B2 respectively. All 42 small subunit 
ribosomal RNA gene sequences were used for 
this phylogenetic study were listed in Table 1. 

Results 

The expected amplicons with sizes of 770 and 
90B base pairs (bp) were observed in all of the 
examined samples (Fig. 1). The sequences of the 
IBS rRNA gene of the T. annulata Iran vaccine 
strain and Iraq (Duhok isolate) were determined 
fmm. the overl^ping flanking sequences of two 
generated PGR fragments. The sequencing of 
the PGR products yielded 1424 and 1413 base 
pair length for two Iran and Iraq samples re- 
spectively. The sequences were then subjected 
to phylogenetic analysis by using the BioE dit 
programme. The sequence identity matrix was 
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also determined and showed 88.2-100 % ho- 
mology between T. annulata 18S rRNA sequence 
of S15 Iran vaccine strain and Thalaia spp. from 
all over the world (Table 1). The sequence of 
Iran vaccine strain showed completely same as 
Iraq sequence (100% identical). But the similar- 
ity of Iran vaccine strain with other T . annulata 
reported 18S rRNA gene sequences from China, 
Spain and Italy determined the 97.9 to 99.9% 
identity. The closest similarity between T. annu- 
lata Iran vaccine strain and otiier Thalffia species 
was belonged to T. lestoquardi Iran vaccine strain 
(AF081135) 99.5% simHarity. There were s^ 
lected 6 outgroup 18S rRNA gene sequences 
from Babesia ovis, B. motasi, B. cabali, B. ivffps, 
Toxoplasma gondii and Isospora suis those exhibited 
the most difference among the comparison 
(88.3-88.7% simHarity). 



1 2 M 




Fig. 1 Gel agarose (2%) electrophoresis of ampli- 
fied 18S rRNA gene sequence of Thalffia annulata. 
Two fragments of 18S rRNA gene sequences were 
amplified for Thalffia annulata. Lane 1; fragment of 
908 bp, lane 2; fragment of 770bp and, M; 100 bp 
D N A ladder as size marker 

Entropy plot 

The alignment was manually edited in BioEdit 
software and truncated to the size of the smal- 
lest sequence (1,413 bp). After resizing the 
aligned 44 sequences, in order to have the cor- 
rect comparison, the designed entropy plot was 
plotted and this design showed the amplified 
1413 bp length of the 18S rRNA gene sequence 
in this study, spans most of the hypervariable 
regions exist along the alignment (Fig. 2). 
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Fig. 2: Plotted entropy shows the 4 major hyper- 
variable regions along the 18S rRNA gene se- 
quence. Two studied T. annulata sequences had 
141 3bp length. The plot shows the fragment of 
Thaleria annulata have an enough length to cover 
the majority of hypervariable regions. The boxes 
are the hypervariable regions along the gene se- 
quence and the arrows demonstrate the region was 
sequenced and ^plied for phylogenetic analysis 

Sequence alignment and phylogenetic 
analysis 

The phylogenetic tree was constructed based 
on the T . annulata Iran SI 5 vaccine strain, Iraq 
isolate, 36 Thalffia spp. and 6 out-group se- 
quences including Babesia spp.. Toxoplasma gon- 
dii and Isospora suis sequences (Table 1). 
Thalffia sequences was divided into six clades 
in the constructed tree, including "Q" clade; 
consists T. equi, "Y" clade; includes T. youngi 
and T. bioomis, "0" clade; contains T. ovis, T. 
capreoli and T. vdifffa, "M" clade; includes T. 
mutans, "B" clade; consists T. buffdi, T. sffgaiti, 
T. sinaisis, T. sqjarata and T, oriaitalis, "A" 
clade; includes T. annulata, T. parva, T. lesto- 
quardi and T. taurotragi, and two more clades 
for outgroup sequences; "Babesia" clade; con- 
sists B. ovis, B. divffgais, B. motasi and B. cabah, 
and "Coccidian" clade; contains Toxoplasma 
gondii and Isospora suis (Table 1 and Fig. 3). 
There are eight T. annulata 18S rRNA gene 
sequences in clade "A", identity percent of T. 
annulata Iran vaccine strain and other T . annu- 
lata in this clade are 99.2% to 100%. 
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Table 1 The sequence identity values between T. 
annulata Iran SI 5 Vaccine Strain and 43 small sub- 
unit of Ribosomal RNA gene sequences belong to 
different Theilffia species and four Babesia spp. as 
well as two outgroup coccidians, Toxoplasma and 
Isospora. Different Theilffia species are categorized in 
six clades (Q, Y, 0, M, B and A) ia addition to two 
outgroup clades for Babesia spp. and Cocddian spp. 



Sequence 


Accession 
number 


Clade* 


Identity 
percent" 


T.equi 


AB515315 




90.1 


T.equi 


AY 534882 






T.eqiii 


r,U 00091)0 B 






T.yoimgi 


AF245279 




91.1 


T.bicomis 


AF499604 


Y 


92.7 


T.velifera 


AF097993 




95.3 


T.capreoli 


AY 726011 




96.4 


T.ovis 


AY 260171 




97.3 


T.ovis 


AY 260172 




97.3 


T.ovis 


EU622911 


- 0 


97.1 


T.ovis 


AY 533144 




97.0 


T.mutans 


AF078815 




94.2 


T.mutans 


FJ213586 


M 


92.9 


T.buffeB 


DQ 287959 




95.7 


TAuffeli 


EF126184 




96.1 


TJbuffeH 


AF236094 




95.9 


TJiuffeH 


DQ 104611 




95.9 


TiuffeH 


Z 15106 


- u 


95.9 


T.bufte]i 


AF236097 


D 


95.9 


T.biifte]i 


FJ426360 




96.6 


T.sergenti 


GU 143088 




92.1 


T.serc(enti 


FI225392 




95.5 


T.sergenti 


EU083803 




95.8 


T. separata 


AY 260175 




95.4 


T.orientalis 


AB520958 




95.9 


T. sinensis 








T.parva 


L02366 




98.5 


T.taurotragi 


L 19082 




97.9 


T. annulata 


EU073963 




99.8 


T. annulata 


EU083799 




99.4 


T. annulata 


EU083800 




99.8 


T. annulata 


EU083801 


- A 


99.9 


T .annulata 


FJ426369 




99.2 


T. annulata 


DQ 287944 




99.6 


T. annulata, 
Iran Vacdne 


HM628581 




ID 


Strain 








T. annulata Iraq 


HM628582 




100.0 


T.lestoquanli 


AF081135 




99.5 


B.ovis 


AY 150058 




85.3 


B.motasi 


AY 533147 


Babesia 


87.0 


B.divergens 


AY 572456 




88.7 


B. cabal! 


AY 534883 






Toxoplasma 
gondii 


L37415 


Coccidian 


84.4 


Isospora suis 


U97523 




83.3 



, ... . 




a 



* The analyzed sequences were grouped in eight major 
clades, based on 18S rRNA gene sequences 
** Nucleotide identities are given in percentage 



Fig. 3: Unrooted phylogenetic tree is constructed 
using Thaleria spp. 18S rRNA gene sequences. 
Thaleria spp. sequences were grouped in six m^or 
clades (detailed description is in the text). Scale bar 
represents nucleotide substitutions per position 

The identity percent between T . annulata Iran 
strain and T. lestoquardi was 99.5% and this 
similarity rate was more than the identity per- 
cent between T . annulata Iran and Italy or Chi- 
na, Xinjiang isolates. 

Discussion 

Based on available literature, this is probably 
the first phylogenetic analysis, molecular char- 
acterization of T. annulata infected cell line of 
Iran vaccine strain, and a field isolate from 
chnically infected cattle in Iraq, Duhok by us- 
ing 18S rRNA gene sequence. 
One of the most important points in Good 
Manufacturing Practice (G MP) is to character- 
ize the local vaccine seed in Razi Vaccine and 
Serum Research Institute, Karaj, Iran. There- 
fore, we decided to classify the T . annulata Iran 
vaccine strain cell line. Hence, we focused on 
to establish the phylogenetic relationships of 
Iran T. annulata vaccine strain with other Thsl- 
aia species using 18S rRNA gene sequences. 
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Basically, there are three steps in all phyloge- 
netic analysis; multiple alignment of the se- 
quences, distance calculation and tree con- 
struction. Using taxa (the outgroup) that are 
known to fall outside of the group of interest 
(the ingnoup) is the way to root tree. In this 
study, we applied Babesia species and two coc- 
cidian parasites to rootThalaia species. 
In addition, Thdleria species are host and vec- 
tor specific (26, 27) but in some countries 
more than one species can infect animals, 
which causes a problem in diagnosis and epi- 
demiology. Although recent molecular studies 
suggest that the genus Thaleria including T. 



annulata, T. parva and T. lestoquardi are very 
similar in features of microscopic characteris- 
tics; these species are phylogenetically distinct 
and can be differentiated by accurate molecu- 
lar techniques (28, 29). Although, Thdlaia spp. 
are tick-transmitted and the parasite is fully 
corresponded to the specific hosts, and this 
association makes a geographical distribution 
of the Thaleria species around the world, but 
the phylogenetic analysis clearly shows close 
relationship of different species, T. annulata, T. 
lestquardi, T. parva and T. taurotragi in con- 
structed tree all are within clade "A" (Fig. 4). 




r-.-i>i'::Tim no 



Fig. 4: Rooted phylogenetic tree of gene sequences of Thaleria annulata S15 Iran vaccine strain, Iraq isolate 
and other Thaleria spp. also outgroup 18S rRNA gene sequences existing in G enBank (Outgroup sequences are 
the least related to the group of taxa that we are studying). The outgroups are considered for rooting a tree. 
Scale bar represents nucleotide substitutions per position 



In this study, two fragments of 18S rRNA 
gene sequences of T. annulata 770bp and 
908bp were amplified for two examined sam- 



ples. Nucleotide sequence identity data dem- 
onstrated that the Iran S15 vaccine strain 
(G enBank Accession Nr^ HM628581) has nu- 
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cleotide homology of 99.9% with China, Nei- 
meng (EU083801), 99.8% with China, Ningxia 
(EU083800), 99.8% with China, Yining 
(EU073963), 99.2% with Italy (FJ426369), 
99.6% with Spain (DQ 287944), 99.4% with 
China, Xinjiang (EU083799) and 100% with 
Iraq, Duhok isolate (HM628582). From this 
perspective, the SI 5 live attenuated vaccine 
strain has a great homology with other T . an- 
nulata isolates; Iraq in neighborhood and Chi- 
na, Spain and Italy far from Iran. But the ef- 
fectiveness of this vaccine for use in these 
countries as a valuable vaccine is an objective 
of future studies. 

The minimum identity belongs to T. annulata 
Italy isolate and maximum similarity was for 
Iraq isolate. The interesting finding is that all 
pathogenic Thalffia species including T. annu- 
lata, T. parva and T. lestoquardi are in "A" clade 
(Table 1). 

Furthermore, T. lestoquardi and T. annulata 
share an extremely high similarity (99.5%) in 
their 18S rRNA gene sequences. This finding 
and the recent occurrence of T . lestoquardi in- 
fection in bovine cells might be a good expla- 
nation for very high closeness of these two 
Thalffia species. 

In the present study, another sequence identity 
analysis was performed by using a smaller size 
of 18S rRNA gene sequence of 793 nucleo- 
tides. The T . annulata field isolate from G oles- 
tan Province of Iran (HM535613) showed the 
100% homology with the T. annulata Iran S15 
vaccine strain, Iraq Duhok and Turkey isolates 
(data not shown). 

Molecular phylogeny was used to classify Thei- 
leria spp. obtained in the Hubei province of 
China; the results clearly determined that the 
Thalffia spp. from ruminants found in Hubei 
belonged to the benign group of Theileria spp. 
(30). 

Three Theileria genotypes were phylogeneti- 
cally analyzed by Nagore et al. in Spain that, 
sharing 96.7%- 97.0% similarity between their 
18S rRNA gene sequences: T. ovis, Theileria sp. 
OTl,andThaleriasp. OT3(29). 



Yin et al. in China suggested that the infective 
Thalaia for small ruminants are comprised of 
two species with the same morphology, vector 
and life cycle. Regarding to the classical tax- 
onomy, they should be one species, but the 
phylogenetic tree showed several Theilaia spp. 
including Chinal and 2 that were isolated 
from sheep and goats (31). 
Nijhof et al. showed the cause of deaths in 
wild ruminants in South Africa that was pre- 
viously attributed to the genus Cytauxzoon, 
but phylogenetic analysis has clarified the Thei- 
kia is the cause of mortality (32). 

In conclusion, according to the data presented 
here, there are high homologies between T. 
annulata Iran 18S rRNA gene sequence strain 
with other T . annulata from above mentioned 
countries, in particular Iraq and Turkey. 
Therefore, if this phylogenetic data correlate 
with immunological response of susceptible 
cattle, it might be a new sight to find an effi- 
cient vaccine to control and prevention of 
Tropical Theileriosis through molecular epi- 
demiological methods. 
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